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Abstract. Neutrino-induced recoil events may constitute a background to direct 
dark matter searches, particularly for those detectors that strive to reach the ton- 
scale and beyond. This paper discusses the expected neutrino-induced background 
spectrum due to several of the most important sources, including solar, atmospheric, 
and diffuse supernova neutrinos. The largest rate arises from produced solar 
neutrinos, providing upwards of 10'^ events per ton-year over all recoil energies 
for the heaviest nuclear targets. However the majority of these *B events are expected 
to be below the recoil threshold of modern detectors. The remaining neutrino sources 
are found to constitute a background to the WIMP-induced recoil rate only if the 
WIMP-nucleon cross section is less than 10""'^^ pb. Finally the sensitivity to diffuse 
supernova neutrino flux for non-electron neutrino flavors is discussed, and projected 
flux limits are compared with existing flux limits. 
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1. Introduction 

Weakly-Interacting Massive Particles (WIMPs) are a leading candidate for the dark 
matter that constitutes ~ 23% of the mass density of the Universe [Ij. A promising 
means to search for WIMPs is through direct interactions in underground detectors, 
whereby the struck nucleus recoils coherently and its kinetic energy is deposited into 
scintillation light. At the present, detectors constrain the spin-independent WIMP- 
nucleon cross section to be S9 x 10~^^ cm^ for WIMPs in the mass range SlOO GeV 
[21, [3] . Future detectors are expected to reach the ton-scale and beyond, providing the 
exciting prospect of covering much of the model parameter space for leading dark matter 
candidates [4j. 

Present direct detection experiments report upper limits on the WIMP mass and 
cross section under the assumption of unknown backgrounds [5j. As the sizes of these 
detectors continue to increase, backgrounds that are beyond the sensitivity limit of 
modern detectors will reveal themselves, necessitating a detailed understanding of these 
background spectra and new techniques for their reduction. Due to the nature of the 
detected events, background subtraction requires understanding the rate and spectral 
shape of any signal that leads to a coherent nuclear recoil event in the energy window 
where a dark matter signal appears. 

In his seminal work over thirty years ago, Freedman [6J pointed out that the 
neutrino- nucleon neutral current interaction leads to a coherence effect, whereby the 
neutrino-nucleus elastic scattering cross section is enhanced and scales approximately 
as the square of the number of neutrons in the nucleus. Coherent recoils are expected if 
the three- momentum imparted to the nucleus, g, obeys the inequality 1, where R is 
the radius of the nucleus. For typical nuclear radii, coherent scattering leads to nuclear 
recoils in the range of a few keV for incoming neutrinos with energies in the range 
~ 1 — 100 MeV. Though neutrino-nucleus coherent scattering remains a fundamental 
prediction of the Standard Model, this process has yet to be observed experimentally 0. 

Solar, atmospheric, and reactor neutrinos may induce observable nuclear recoil 
events [8j. Monroe and Fisher ^ have recently determined the coherent scattering 
rates in future dark matter detectors due to solar, atmospheric, and geo-neutrinos, in 
particular identifying ®B solar neutrinos as the source with the largest event rate. For 
detectors at the ~ ton-yr exposure scale, ®B neutrino-induced recoils are significant for 
recoil energies S3 keV for heavy targets such as Xe and S30 keV for lighter targets 
such as Flourine and Carbon. Assuming a ~ 100 GeV WIMP, neutrino-induced 
recoils constitute a background to the WIMP recoil spectrum if the WIMP-nucleon 
cross section is £ 10~^° pb [9llT0]. 

The goal of this paper is to provide a survey of the expected event rate recoil 
spectra from several important astrophysical neutrino sources. Solar, atmospheric, and 
diffuse supernova neutrinos are studied; geo-neutrinos are not considered since their 
signal is dwarfed by the solar neutrino signal over a similar energy range |9]. In addition 
to providing an updated calculation of the solar and atmospheric neutrino fluxes for 
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a variety of often-discussed nuclear targets, the first predictions for the recoil energy 
spectrum of the Diffuse Supernova Neutrino Background (DSNB) are provided. If the 
neutrinos can be fit, or if the dominant component to their recoil spectrum lies below 
the analysis energy threshold, astrophysical neutrino backgrounds will only contaminate 
a WIMP-nucleon recoil signal if the WIMP-nucleon cross section is ^ 10~^^ pb. It is 
additionally shown that, with current exposures, detectors should have sensitivity to 
reduce the current upper limit on the and v^. component of DSNB flux by more than 
an order of magnitude. 



2. Neutrino-Induced Recoil Spectra 

2.1. Neutrino-Nucleus Coherent Scattering 

The neutrino-nucleus scattering cross section per recoil kinetic energy T is [H] 

Here = — (1 — 4sin^0^)Z is the weak nuclear charge, where is the number 
of neutrons and Z is the number of protons. The mass of the nucleus is M = AM^v, 
where A = N + Z the mass number and M^r = 931 MeV. From kinematics the recoil 
energy is related to the neutrino three-momentum via = 2MT, and the maximum 
three-momentum transfer is a function of the neutrino energy, q^^x — The range 

of recoil energy for the struck nucleus varies between and Tmax, where 

Taking the limit M S> ^E^,, which will always be valid for the results presented here. 



the minimum neutrino energy for a fixed recoil kinetic energy T is E,^ = yMT/2. 

The form factor, F{Q'^y, is a function of the 4-momentum, Q, and measures the 
distribution of nucleons within the nucleus and thus the deviations from coherence at 
high recoil energies, defined by momentum transfers qR ^1. The form factor at zero 
momentum transfer scales as -F(O)^ —>■ 1 for qR -C 1. To approximate the effects of 
the form factor, the analysis in this paper uses the Helm form factor as parametrized in 
Lewin and Smith [12] . As in the case of WIMP-nucleon recoils, the form factor has the 
primary effect of reducing the cross section for large momentum transfers. 



2.2. Neutrino Fluxes 

Since the coherent scattering process is "flavor blind," mixing angles which connect 
neutrino weak flavor eigenstates to mass eigenstates need not be considered in flux 
calculations. For all of the sources, the energy spectrum for the weak flavor eigenstates 
at the production sites are used. 

Solar Neutrinos - The calculations in this paper consider two sources of solar 
neutrinos, so-called '^B and hep neutrinos. These neutrinos come from the reactions 
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®B — i>'^Be* + e+ + z/e and ^He +p — ^''He +e+ + z/g, which occur in 2 x 10~^% and 
2 X 10~^% of the terminations of the solar pp chain, respectively fT3]. The remaining 
sources of solar neutrinos induce recoils with energies and rates below that of the 
induced spectrum. For example, neutrinos produced directly in the pp reaction only 
induce recoils up to a maximum kinetic energy of ~ 0.03 keV [9|. 

Atmospheric Neutrinos- Cosmic ray collisions in the atmosphere produce a primary 
flux of i^e, t'e, i^/i, and |i4j, with a spectrum that extends up to energies of ~ 100 
GeV. At these high energies, the loss of coherence in the scatterings must be accounted 
for. On the one hand, the coherent scattering cross section scales as the neutrino energy 
squared, favoring high neutrino energies. On the other hand, the form factor suppresses 
scatterings for high momentum transfers, which largely result from the highest energy 
neutrinos in the flux distribution. Convolving both of these factors means that detectors 
are mainly sensitive to the atmospheric neutrino flux S 100 MeV. A more exact upper 
cut-off depends on the detector target; the larger the mass number of the nucleus, the 
more strongly the event rate is suppressed by the form factor at the highest detectable 
recoil energies. For the atmospheric neutrino flux, the calculations in this paper use the 
published Fluka results for the Kamioka location, which have been determined down 
to neutrino energies ~ 10 MeV [15]. At energies of interest for this calculation, £100 
MeV, there is a ~ 20% uncertainty in the normalization of the flux, primarily due to 
geo-magnetic effects. 

Diffuse Supernova Neutrinos- Core-collapse supernova produce a burst of ~ 10^^ 
neutrinos of all flavors per explosion; the past history of all of these explosions produce 
a Diffuse Supenova Neutrino Background (DSNB). The DSNB flux is a convolution of 
the core-collapse supernova rate as a function of redshift with the neutrino spectrum 
per supernova. The core-collapse rate is derived from the star-formation rate and stellar 
initial mass function; the calculations in the paper follow the recent results of Horiuchi et 
al. [H] . The neutrino spectrum of a core-collapse supernova is believed to be similar to a 
Fermi-Dirac spectrum, with temperatures in the range 3-8 MeV [T7|. The calculations 
in this paper assume the following temperatures for each neutrino flavor: Ty^ = 3 MeV, 
Tp^ = 5 MeV, and T^^ = 8 MeV. Here T^^ represent the four flavors u^, u^, u^-, and Ur- 
Due to the ~ E'^ scaling of the total cross-section in Eq. [T] (integrated over all recoil 
energies), the flavors with the largest temperature dominate the event rate. 

There are strong upper limits on the component of the DSNB flux from Super- 
Kamiokande, < 1.2 cm^^ s^^ [IE]. This flux limit constrains the high energy tail 
of the DSNB flux, specifically neutrino energies ^ 18 MeV. More recently, the SNO 
collaboration has placed an upper limit on the z/g component of the flux of ~ 70 cm~^ 
s~^ for neutrinos in the energy range of 22.9 — 36.9 MeV [19]. Though weaker than 
the z/g limit, the SNO z/g limit does probe the parameter space for exotic core-collapse 
models [20]. An indirect bound on the z/g flux that rivals the present limit on the z/g flux 
follows from oscillation conditions [21], and stronger direct bounds on the z/g flux may be 
realized in future Argon TPC detectors [22]. The strongest bound on flavors comes 
from the Mt. Blanc experiment [23], which quotes an upper limit on the flux of ~ 10^ 
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cm~^ for neutrino energies 20 MeV. This flux limit was set based on the neutral 
current interaction with ^^C. It has additionally been proposed that Super-Kamiokande 
can set an improved limit on the Vx flux [2^. Below the prospects for improving on 
these flux limits are discussed. 

In addition to the DSNB, detectors sensitive to low energy nuclear recoils may be 
used to detect a supernova burst in the Galaxy [25]. For exposures on the scale of 
~ ton-yr, there are expected to be of order 10 events, summed over all flavors for a 
supernova at 10 kpc [26] . The detection of a supernova burst via a recoiling nucleus 
would provide a unique window into supernova physics, complementing the detection of 
proton recoils [27] as the best means to study the component of the neutrino flux. 

Figure [1] shows the neutrino flux spectra for solar, atmospheric, and DSNB sources. 
In terms of integrated flux, the ®B flux is by far dominant, nearly three orders of 
magnitude larger than the next closest hep rate. However, as the flux is confined to 
energies ^ 16 MeV, these will produce lower energy recoils in comparison to the higher 
energy DSNB and atmospheric components. A convolution between the flux and cross 
section is required to obtain the recoil rate as a function of energy; the goal of the 
following section is to determine these respective rates. 




Neutrino Energy [MeV] 



Figure 1. Neutrino fluxes for astrophysical sources that constitute tlie dominant 
backgrounds to WIMP recoil signals. From top to bottom on the left, the fluxes 
are solar *B (solid, red) and hep (dotted, blue), DSNB (long-dash, black) with 
temperatures of 3 5, and 8 MeV (corresponding to I'e, and v^, respectively, where 
Vx symbolizes muon and tau neutrinos and the respective antiparticles). The short- 
dashed (green) curves at the highest energies are the atmospheric neutrino fluxes, 
plotted down to the lowest energy in the calculation of Ref. 15 . From bottom to top 
on the right, the fluxes are for v,,, v^, then i/^, D^. 
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Figure 2. Event rate per recoil kinetic energy for six target nuclei. For both the 
diffuse supernova and atmospheric event rates, the sum of all contributing neutrino 
flavors are shown. 
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Figure 3. Number of events above a threshold recoil kinetic energy for six target 
nuclei. For both the diffuse supernova and atmospheric event rates, the sum of all 
contributing neutrino flavors are shown. 
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2.3. Event Rates 



The differential event rate at a fixed recoil kinetic energy is 
dR{T) _ dN da{E^,T) 
dT JEniin dEi, dT 



(3) 



where dN/dE^ is the neutrino flux spectrum and -Emin = \JMT/2 is the minimum 
neutrino energy for a given recoil energy as dictated by the kinematics. Figure [2] 
shows the event rate recoil spectrum for six different nuclear targets. For all targets 
the naturally-occurring abundances are assumed. For both the diffuse supernova and 
atmospheric event rates, the sum of all contributing neutrino flavors are shown. In 
particular for the DSNB, an 8 MeV spectrum from Figure [T] is multiplied by four to 
account for the production spectrum of the four i>x flavors. Due to their relatively hard 
spectra, the flavors are seen to dominate the event rate, particular at high recoil 
energies; there is only about a ~ 10% increase by including the 3 and 5 MeV spectra 
at the lowest recoil energies. Each of these curves are the true, infinite resolution 
spectra, i.e. they do not account for the expected finite energy resolution of detectors. 
A detailed convolution with a resolution function will depend on the nuclear target and 
the particular experimental environment. 

Figure [3] shows the number of events above a given recoil energy for the same six 
nuclear targets. Most of the events are confined to low recoil energies, for example 
for the case of Xe there are a total of ~ 10'^ events over all energies, but only ~ 1 
event per ton-yr above a threshold of 3 keV. Future Xe detectors are expected to have 
thresholds in the area of ~ 5 keV; as is seen dropping the threshold below this energy 
will lead to a significantly increased signal. 

As an additional note, the analysis above just accounts for neutrino-nucleus 
coherent scattering. In principle it would also be possible to detect these same fluxes 
via neutrino-electron elastic scatterings [8]. For this channel the largest rate is to due 
the solar pp reaction. For example, from pp scatterings on Xe a flat spectrum of electron 
recoil events is expected at ~ 0.1 events per ton-yr with energies up to ~ 600 keV. 

3. Implications for WIMP-Nucleon Cross Section Constraints 

In the absence of backgrounds the expected upper limit on the WIMP-cross section 
simply scales linearly with the detector. For example a ten times greater exposure 
will imply a ten times stronger upper limit on the cross section. In the presence of 
backgrounds, however, the projected limits on the cross section must be modified. 
Dodelson [28] has provided a simple formalism for estimating the upper limit on the 
WIMP-nucleon cross section, given a measured background rate and a fiducial detector 
volume. In this formalism, the probability of observing a total of events, given a 
WIMP-nucleon cross section, a, is 
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Here the mean, /x, is a sum of the WIMP-induced signal and neutrino-induced 
background events, and depends on the model parameters such as a and the expected 
background spectrum. The mean number of background events is iV^ and the error on 
the background is a^. In Eq. Hlthe effect of binning the data in energy is neglected, and 
as a result the estimates quoted are likely to be conservative (binning the data would 
likely increase the sensitivity). 

From the likelihood in Eq. HJ Dodelson [28] defines the Background Penalty Factor 
{EPF) as 

upper limit 

upper limit '^current ^ , r> d / c \ 

f^future = >< BPF, (5) 

where A is the ratio of the future detector exposure to the current exposure. The EPF 
then essentially measures, for a given background rate, the projected upper limit on the 
cross section for a future detector. It is calculated by integrating the normalized version 
of Eq. m up to the desired confidence level, and the projected cross section upper limit 
is then obtained from Eq. [51 Here the projected upper limits are quoted at 95% c.l. 

The results presented here assume that the mean background event rates are 
obtained from the results in the figures above. The error on the background is taken 
to be (Jh = 0.2iVfe, motivated by the uncertainty on the solar neutrino fluxes [29]. The 
results presented here are found to be independent of similar reasonable assumptions for 
the error on the background. As far as the actual measured backgrounds are concerned, 
recent results of XENONIO [2] and CDMS [3] give background rates per unit volume of 
(7 events)/(0.44 kg-yr) and (0.6 events)/(0.33 kg-yr), respectively. These backgrounds 
are not accounted for in this present analysis, amounting to the assumption that they 
can be reduced and understood in future experimental analysis. The calculations above 
show that, given the current exposures, the total background rate from all neutrino 
sources is ~ 10~^ events above CDMS and XENONIO thresholds, so it is unlikely that 
these background events are neutrino-induced. 

Accounting for the neutrino backgrounds. Figure H] shows the upper limit on the 
cross section attainable as a function of detector exposure. Two target nuclei are 
considered, Xe and Ge. To set the absolute scale, the current WIMP-nucleon cross 
section upper limits from these experiments are ~ 4 x 10~^^ cm^ (at ~ 30 GeV) and 
~ 5 X lO""^*^ cm^ (at ~ 60 GeV), respectively [21 [3j. Each target shows the results 
for three different energy thresholds. The salient point is that, as the threshold is 
increased, the dominant ®B background is reduced and the projected upper limit on 
the cross section will become much more stringent. In fact, cutting out the signal 
entirely (which amounts to adopting the 5 keV threshold for Xe and the 7 keV threshold 
for Ge), the projected upper limit will continue to decrease to below ~ 10~^^ pb, the 
lower regime of the supersymmetric model parameter space [4]. 

It should be noted that, while effective at eliminating backgrounds, setting the 
threshold at too high of an energy may have the adverse effect of cutting out the sought 
after signal events from WIMP-induced recoils. For example, for a 100 GeV WIMP, 
~ 10% of the recoil events occur in energy bins below 3 keV. Similarly for a 50 GeV 
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WIMP, ~ 15% of the events occur in energy bins below 3 keV. So setting the threshold 
to cut the entire rate may not be the optimal strategy, depending on the underlying 
WIMP parameters. In this sense, determining the WIMP signal robustly will require a 
multi-component fit to both the WIMP and neutrino-induced recoil spectra. 
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Figure 4. The upper limit on the WIMP-nucleon cross section, relative to the current 
upper limit, as a function of the exposure, relative to the current exposure. Two 
target nuclei are shown, Xe and Ge. The labels in each case denote the lower energy 
threshold assumed for the experiment. It is assumed that there are no non- neutrino 
sources contributing to the background. 



4. Upper DSNB flux limit on non-electron neutrino flavors 

The preceding section has discussed the extent to which the neutrino backgrounds can 
degrade the upper limits on the WIMP-nucleon cross section. It is also of course 
interesting to explore the neutrino fluxes in the context of a signal themselves. The 
focus here is on the DSNB flux, and more specifically to provide an estimate for the 
expected flux limit that current detectors can set on the flux component. Similar 
to the preceding section, the projections focus on Xe and Ge targets, though they can 
easily be extended to other targets with similar results. 

From Figure [3l the expected DSNB event rate is ~ 0.05 events per ton-yr. As 
discussed above, this event rate is dominated by the four flavors, again not accounting 
for neutrino mixing. It follows that the detection of this flux requires an exposure of 
the order 100 ton-yr. To set the current scale, given that the current exposure for Xe is 
0.44 kg-yr, the expected number of DSNB events is a paltry ~ 2 x 10~^. 

Nonetheless, from these predicted event rates an estimate of the flux upper 
limit that present detectors can set may be obtained by simply increasing the flux 
normalization so that 0{1) event would be expected in detectors with the present 
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exposures. The implied increase in the flux normalization is ~ 1/(2 x 10~^) = 5 x 10^. 
Given that the current best estimates for the DSNB flux are ~ 10 cm~^ s~^ for non- 
electron flavors (again not accounting for neutrino mixing effects), it is projected that, 
with current exposures, the upper flux limit will be reduced by more than an order 
of magnitude, so that the projected upper flux limit is 

^limit <^Q6^^~2g-l_ (g^ 

This limit is valid for all energies, so in this way it is better than an order of magnitude 
greater than the Mt. Blanc flux limit, which is only valid for neutrino energies ^20 
MeV. 

The flux upper limit in Eq.[6]may be viewed as optimistic, because no contaminating 
backgrounds in the energy window of the DSNB signal have been accounted for. 
Accounting for the current backgrounds may somewhat degrade this limit, though as 
in the preceding section it is assumed that these backgrounds can be eliminated in 
future runs of the experiment. Going further and pushing the flux sensitivity even 
lower requires subtracting the solar and atmospheric backgrounds. In fact, in order to 
determine how much better the flux limit in Eq. [6] will get as a function of detector 
exposure, it is possible to again appeal to the EPF formalism from the previous section. 
The only difference is that replaces the cross section as the quantity to constrain. 

Under the assumptions above the projected upper limits on the DSNB flux are 
obtained simply by reading off the limits from the curves in Fig. HI Since the DSNB 
constitutes a small fraction of the total backgrounds in Fig. HJ it is found that calculating 
the flux limit in this manner provides an excellent approximation to the case in which the 
curves in Fig.Hlare calculated without the DSNB fiux included. To set the normalization, 
the current flux limit is taken to be the projected value of (/>[l™^* from Eq. [61 For a 
keV threshold using an Xe target, from Fig. H] the upper flux limit is projected to be 
reduced by about an order of magnitude further to ^ 10^ cm~^ s~^ for a two order of 
magnitude increase in exposure. This scaling is again primarily because of the solar 
neutrino background. Eliminating the *B background with an energy threshold cut. 
Fig. m shows that the flux sensitivity will increase linearly with exposure down to the 
projected DSNB flux window of ~ 10 cm"^ s^^ for exposures ~ 10^ times the current 
exposures. 

Finally, it should be noted again that as in the case for the projected cross section 
limits, the DSNB flux estimates have not accounted for the detailed spectral shape of 
the signal or backgrounds. Only the energy threshold has been varied to eliminate the 
dominant ^B background. Future detectors, with orders of magnitude more exposure 
than at present should be able to improve on the upper limits provided here by using 
the difference between the spectral shapes of the signal and background. 

5. Discussion and Conclusion 

This paper has discussed the neutrino-induced recoil spectrum at dark matter detectors 
due to solar, atmospheric, and diffuse supernova neutrinos. It is shown that ton-scale 
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detectors can expect ~ 100 '^B neutrino events per ton year over all recoil energies. The 
dominant signal can be excluded using either spectral information or by an energy 
threshold cut. It should also be noted that the *B signal could be cleanly identified by 
a detector with directional sensitivity [30j . 

The effect of these neutrino-incuded recoils on the extraction of an upper limit on 
the WIMP-nucleon cross section has been examined. If the ^B signal can be excluded, 
and for the particular case of Xe and Ge targets, the upper limit on the cross section is 
projected to scale linearly with the detector exposure. Non- *B backgrounds only impact 
the WIMP-nucleon cross section upper limit if the cross section is below ~ 10^^^ pb. 
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